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ABSTRACT 



Experiments have been conducted in order to determine the 
burning rates of single droplets of liquid fuels in oxidizing atmos- 
pheres. N-heptane and absolute ethyl alcohol were tested in nitrogen- 
oxygen mixtures for oxygen weight fractions between 0. 10 and 0.90. 
Measurements were made also for benzene and toluene; for oxygen 
weight fractions greater than for air, a carbon residue formed during 
the burning of these fuels, which made it impossible to obtain accurate 
data. For all nitrogen- oxygen mixtures it was found that the mass 
burning rate of single droplets of fuel is proportional to the first 
power of the drop diameter. 

The theoretical analysis on the burning of single droplets of 
fuels in an oxidizing atmosphere by Goldsmith and Penner was used to 
obtain values for the evaporation constants of the fuels tested. The 
experimentally determined values of the evaporation constants were 
found to be in good agreement with the calculated values for all fuels 
tested which did not form carbon. 
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I. INTRODUCTION’ 



During the past few years there has been increased effort to 
determine the mechanism of combustion of individual droplets of 
fuel. These basic studies were initiated in order to aid in the design 
of efficient spray injection type burners for heterogeneous combustion. 
Heterogeneous combustion is of importance in such widely different 
applications as stationary boilers, Diesel engines, gas turbines, and 
rocket motors. 

iost previous experimental work on the combustion of liquid 

droplets of fuel has been limited to the use of air as the oxidizing 

( 12 ) 

medium. For example, G.A.E. Godsave ' ' ' determined experi- 
mentally the burning rates of hydro-carbon fuels such as benzene, 
ethyl alcohol, and normal heptane at atmospheric pressure, using 
air as the oxidizing medium. Godsave estimated also the diameter of 
the flame front surrounding the liquid drop for a number of fuels. 

The most important conclusion derived from this experimental work 
was that the mass evaporation rate of the drop is proportional to the 
first power of the drop diameter. 

{ 3 } 

A. R. Hall and J. Diedericksen' ' have carried out experimental 
work on the combustion of drops of liquid fuels. They verified the 
result that the mass rate of burning is proportional to the first power 
of the drop diameter. Experiments were conducted also to determine 
the effect of pressure on the combustion of droplets of fuel. Data were 
obtained for the burning rate at pressures up to twenty atmospheres. 
The results showed that the mass rate of burning of fuel droplets is 
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roughly proportional to the one-fourth power of the pressure. 

In conjunction with the experimental studies, there has plso been 

an effort made towards a theoretical understanding of the mechanism 

of the combustion of single dror>3 of fuel in an oxidizing atmosphere. 

( 12 ) 

Godsave' ' ’ obtained a useful theoretical interpretation of hi3 experi- 
mental results on the assumption that the chemical reaction rates dc 
not control the rate of burning. If this is the case, the problem is 
reduced to a study of heat and mass transfer between the flame front, 
the fuel droplet, and the surrounding atmosphere. The flame front is 
formed at some distance from the drop surface. The oxidizer is 
delivered from the surrounding atmosphere to the region of active 
combustion by convection and diffusion; the fuel evaporates and 
diffuses, without chemicai change, to the reaction front, which is 
assumed to be a spherical shell surrounding the droplet. The location 

of the reaction front is determined empirically. Utilizing the physical 

( 1 2 ) 

model described above,' Godsave' ’ ' was able to derive an explicit 

expression for the burning rate of liquid fuel droplets which contained, 

however, two adjustable parameters, viz. , the temperature of the 

flame front and the combustion radius. Godsave showed that the use 

of reasonable values for these parameters was consistent with the 

observed mass burning rates. 

(4) 

Goldsmith and Penaer ' 1 extended Godsave's work by using the 
hypothesis that the position of the flame front is established in such a 
way that the ratio of fuel to oxidizer corresponds to stoichiometric 
proportions. In this way, they were able to obtain explicit relations 
for the radius of the flame front, the flame temperature, and the mass 
burning rate, which contained no adjustable parameters. As the result 
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of an improved method of formulating the problem, Goldsmith and 
Penner were able to solve the equations without introducing invalid 
approximations concerning the physico-chemical parameters. 

There are obvious deficiencies in the assumed model for a 
burning droplet, which is shown schematically in Fig. 1. The most 
serious deficiencies are: 

1. Contrary to the assumption that a spherical reaction surface 
exists, photographs of burning droplets indicate that the luminous 
region is of the form shown in Fig. 2. Godsave measured the diameter 
of the sphere corresponding to the dotted circle in Fig. 2. The 
existing flame shape is the result of natural convection currents during 
burning, which have been assumed to be negligible in the theoretical 
analysis. 

2. The assumption that no chemical reactions occur in non- 

stoichiometric mixtures is obviously invalid, although a similar 

hypothesis has been used successfully in the description of gaseous 
(5) 

diffusion flames ' . 

3. The calculations of mass burning rate neglect the occurrence 
of dissociation at the highest temperatures, followed by reassociation 
at lower temperatures. A more complete study, without an assumed 
reaction surface, will probably not lead to difficulties of this sort. 

Since previous experimental work on the burning of single drops 
of liquid fuels has been confined to air as the oxidizing medium, the 
comparison between theory and experiment is not as stringent a test 
of the theory as might be desired. With this in mind the author has 
performed experiments on a limited number of fuels in order to 

I 
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determine the burning rate o£ single droplets in nitrogen-oxygen 

mixtures containing varying amounts of oxygen. The evaporation 

constants (burning rates) determined experimentally are compared 

in Section IV of this report with theoretical values for the evaporation 

( 4 ) 

constant obtained from the equations of Goldsmith and Penner' ' . 
Although unbelievable high values are calculated for the flame tempera- 
ture, the observed burning rates are again in accord with the theoreti- 
cal results, thereby suggesting the occurrence of compensating errors 
in the analysis. 
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H. rHEOR LTICi'i L CALCULATIONS 
A. Assumptions 

In order to present a clearer picture of the physical model upon 
which the theoretical analysis of Goldsmith and Penner^ is based, 
the important assumptions are tabulated below. 

1. The droplets are spherical. 

2. Convection effects may be neglected. 

3. The flame front surrounding the drop is represented by a 
spherical surface concentric with the drop. All reactions 
take place instantaneously at this surface, at which the 
delivery rates of fuel to oxidizer are in stoichiometric 
proportions. 

4. Steady-state conditions are assumed for fixed droplet 
sizes. This restriction greatly facilitates the mathematical 
treatment. It is reasonable to assume that the solution 
obtained for a fixed size applies to a drop decreasing in 
size when it reaches the radius used in the steady- state 
solution. 

5. The effect of heat transfer by radiation is neglected. 

6. .‘lean values are used, when appropriate, for the 
physical properties. 

7. The temperature of the liquid drop is assumed to be uniform 
and equal to the boiling temperature. 

S. The pressure is assumed to be uniform throughout the 
system. 
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A schematic diagram of an evaporating and burning fuel droplet 
in an oxidizing atmosphere is shown in Fig. 1. The radius of the 
liquid drop is and its temperature is the normal boiling point Tj. . 
The radial distance of the combustion surface from the center of the 
liquid droplet is and its temperature is T . The oxygen-inert 
gas mixture at a large distance from the combustion surface is at the 
temperature T ( . 

3. Equations Necessary to Determine the Evaporation Constant 

(4) 

The equations derived by Goldsmith and Penner for determin- 
ing the flame temperature T , the combustion radius • at *d the 

mass burning rate are summarized below: 



T c 
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( 2 ) 
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The symbols used in Equations (1) through (4) are defined in the 



Table of Nomenclature. 



The linear functional relation between and £ is 

customarily used to obtain the following expression for the variation 

( i 2 4) 

oi’ droplet diameter with time: * 1 

d 2 = do " K't (5) 

Here d is the droplet diameter at time t, d^is the initial droplet 
diameter, and the evaporation constant K' is defined by the following 
relation: 

K ' = 2 liip/ir Vify (6) 

Since K' is independent of YJ , it is a convenient parameter for 
comparing burning rates of different fuels for arbitrary droplet sizes. 

The procedure for calculating the various quantities involves 
the following steps: 

(a) F or suitably chosen values of the physico-chemical 
parameters a value of T is obtained from Eq, (1). 

(b) The computational parameter § is next obtained from 
Eq. (3). 

(c) The quantity *c/v< is obtained from Eq. (2). 

(d) rrip is determined from Eq. (4). 

(e) Finally the evaporation constant K' is calculated from 

Eq. (6). 

Listed in Table 1 are appropriate values of the physico-chemical 
parameters used to obtain theoretical values for the evaporation 



constant. 



TABLE 1. PHYSICO-CHEMICAL PARAMETERS 
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C. Results 

The theoretical calculations for the determination of the combustion 
and temperature (T c ), the radius of the reaction zone ( IQ, ), and the 
evaporation constant (K ') were curried out by Goldsmith and Penner for 
combustion in air (for which the weight fraction of oxygen, has 

the value 0.23). Similar calculations have been performed by the 
author for oxygen weight fractions Y_ of 0. 10, 0. 37, 0. 50, 0. 70, 
and 0. 90. The results o£ these calculations are listed in Table 2, 
together with the known results for air. 

Reference to Table 2 indicates the expected trend of the data, name- 
ly, the calculated rate of burning increases with an increase in the con- 
centration of oxygen. Furthermore, with a decreased concentration of 
diluent inert gases, the combustion surface moves closer to the liquid 
drop. For a weight fraction of oxygen smaller than that for air, the 
opposite behavior is noted, i. e. , for Y^ Q = 0. 10, the combustion sur- 
face moves farther from the drop and the rate of burning ’is decreased. 

The values obtained for the temperature of the combustion zone, 

T fc , are, of course, fictitious. As was stated in the introduction, the 

high flame temperatures are a direct consequence of deficiencies in the 

( 4 ) 

assumed physical model. Goldsmith and Penner'*' noted that proper 
allowance for dissociation and reassociation of the reaction products 
would greatly decrease the value of T c but would probably not change 
materially the temperature gradient near the droplet surface, which 
determines the burning rate. The success of the theory must be as- 
cribed to this type of compensation, at least until a satisfactory theo - 
retical model can be analyzed quantitatively. 
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III. EXP ESI. TV. 1, VOR;< 

A. Apparatus 

The apparatus used for photographing single burning drops 
suspended on a silica filament is shown in Fig. 3. The object of the 
Investigations was to record the decrease in drop size with time during 
combustion. The combustion tank shown in Fig. 3 was utilised so that 
the composition of the oxidizing medium surrounding the fuel droplet 
could be varied. The tank is cylindrical, 42 inches long, with an inside 
diameter of 6-1/2 inches and an outside diameter of 8 inches. It can 
be pressurized to approximately 500 psi and can be heated to a tempera- 
ture of about 1000°F. There are a total of ten windows, 2 inches in 
diameter, that is, five windows on each side diametrically opposed. 

The tank was constructed by M. Goldsmith for experiments on falling 
drops, combustion under pressure, and combustion at elevated 
temperatures. 

A fuel injection system is built into the head of the combustion 
tank. This system is provided with a fuel reservoir, a fuel inlet valve 
for filling the reservoir, and with means for pressurizing the reservoir. 
The fuel system is connected through a ball check valve to a hollow 
brass needle which extends approximately S inches down into the tank. 

A fine silica filament was cemented to the brass needle. This 
filament was of such length that its tip extended vertically down so as to 
be in the line of sight of the top windows of the tank. It was necessary 
to thicken the end of the filament in order to retain the liquid drop on the 
filament. The diameter of the thickened end of the filament was approxi- 
mately 0. 30 mm. 
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The drops were suspenJed on the filament by forcing fuel through 
the brass needle and onto the silica fila nent. The suspended drops were 
found to be between 1. 5 and l. 8 rum in diameter. s can be seen fro.. 
Fig. 4, the suspended drops were reasonably spherical during the 
major part of their life. 

The suspended drops were ignited by means of an electric spark, 

A six-volt storage battery was connected through a spark coil and an 
on-off switch to a spark plug screwed into the head of the combustion 
tank. A welding rod was used to lengthen the center pole of the spark 
plug to a length commensurate with the tip of the filament. Another rod 
was extended down from the head forming a spark gap approximately 
2 cm in length, across the filament tip. This method of ignition was 
found to be satisfactory for oxygen concentrations greater than about 
20 percent in oxygen-nitrogen mixtures. At Y n = 0. 10 it was im- 
possible to ignite a drop for the four fuels tested. 

An electrically driven Beil and Howell ,, Eyemo M 32 mxri movie 
camera was used to photograph the burning drops. The drops were 
photographed in silhouette by providing strong background illumination, 
and using a specially selected lens in order to obtain as large an image 
as possible of the drop on each frame of the film. The lens used was a 
six inch F/4. 5 Cooke Teleking Anastigmat connected to the camera 
through a 10 inch adapter tube. A lens setting of £16 was determined to 
be the most satisfactory setting. The background illumination was 
provided by a 100 watt light bulb shining directly into the tank window 
opposite that into which the camera was focused. Kodak Plus X 32 mr.i 



film was used. 
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The drops were photographed at a camera setting of 24 frames 
per second. The actual camera speed at this setting was checked by 
photographing a calibrated stop watch. The camera timing v/as found 
to be 24.69 frames per second (1 frame per 0.04065 seconds). 

A 3/32 inch ball bearing was photographed at the beginning and 
end of each 100 foot roll of film used. This calibration was carried out 
under the same camera focusing conditions as for the burning drops, 
thereby providing an accurate reference measurement for determining 
the magnification of the photographs and, consequently, the actual 
size of the drops. 

The size of the burning drops was determined by using a technique 
very similar to that described by Godsave^. The film was measured 
with the aid of a 32 mm microfilm reader (to produce further magnifi- 
cation) and a steel scale graduated to l/lOQth of an inch. Two measure- 
ments were made on each frame, namely, the two perpendicular 
diameters inclined at 45° to the major and minor axes in the plane of 
observation. The mean of these two measurements was recorded as 
the "effective diameter" of the drop. 1£ the major and minor axes do 
not differ greatly, as was the case in our tests, then it is easily shown 
that the volume of a sphere with the measured effective diameter is not 
greatly different from that of the prolate spheroid which actually 
corresponds to the shape e£ the drop. 

In most cases, measurements of each drop were taken over a 
range of size extending from ignition down to approximately one half 
or the original drop diameter. T or burning in air, every fifth frame of 
the film on each drop was measured. As the burning rate inure* sed, 
with increase in oxygen concentration, every two to four fra* .cc of the 
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film for any given drop were measured. 

In order to vary the weight fraction of the oxygen (Y, ) in the 

atmosphere inside the combustion tank, a bottle of compressed oxygen 
was connected through a regulator valve to the tank. Then, for example, 
for Y ( „ q * 0,50, the combustion tank was pressurized to a gage 
pressure of 7. 1 psi. The tank was kept at the required pressure for an 
extended time in order to allow complete mixing of the oxygen with the 
original air in the tank. The pressure in the tank was slowly decreased 
to the operating pressure, a drop of fuel suspended on the filament, the 
camera started, and the drop ignited. After each run the combustion 
tank was flushed with air before repeating the above process for another 
run. All combustion experiments were carried out at atmospheric 
pressure. It should be noted that for Yq q = 0. 10 the combustion tank 
was pressurized with nitrogen to a gage pressure of 19.9 psi before 
bleeding to atmospheric pressure. However, ignition at this reduced 
oxygen concentration wa3 impossible and no experimental data were 
obtained, 

B. Experimental Results 

A systematic study was made to determine the evaporation 

constant, X', for different oxygen concentrations and for the following 

fuels: n-heptane, absolute ethyl alcohol, benzene, and toluene. The 

experimental results gave directly the droplet diameter, d, as a 

2 

function of the time, t. The plots of d against t were found to be 

linear. Typical experimental plots for n-heptane and ethyl alcohol for 

Y,, = 0. 50, and for benzene and toluene at = 0.23 are shown in 

O, o O.o 
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Figs. 5 to 8, respectively. Thus the linear functional relation between 

2 2 

drop diameter d and time t, d = d - K't, has been verified again. 

Here d Q is the diameter at time t = o (corresponding to the beginning 

of steady burning) and K' is the evaporation constant. The value of K 1 

2 

is determined directly from the slope of plots of d vs t and has the 
dimensions cm^ sec *. 

The experimentally determined values of K' for the four fuels 
tested at different oxygen weight fractions are listed in Table 3. The 
results of five separate experiments were examined for each fuel at 
each value of Yq q . The tabulated values of the evaporation constant 
are the mean values for each set of five runs. The reproducibility in 
each set of runs was good; the maximum spread of experimental data 
was found to be less than t * percent. 

With the lighting used for photography during the experiments, it 
was impossible to obtain pictures of the flames suitable for measurement 
of flame surface diameter. For n-heptane and ethyl alcohol the luminous 
flame surface was not sufficiently well defined to permit conclusions 
regarding flame shape and flame size. For benzene and toluene burning 
in air, the luminous region could be seen but was not sufficiently clear 
to yield accurate measurement. However, it was possible to determine 
qualitatively that for the latter cases the apparent values of Tc/fo are 
considerably smaller than the theoretical values listed in Table 2. This 
noticeable difference between the observed and calculated values for 

(4) 

Ic/ICe has been considered in the paper by Goldsmith and Penner' , 

and should be ascribed to deficiencies in the theoretical treatment, which 
does not make proper allowance for the influence of convection currents 
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TABLE 3. EXPERIMENTAL VALUES OF K’ (cm 2 sec" 1 ) FOR SEVERAL 
FUELS AND WEIGHT FRACTIONS OF OXYGEN 



Fuel 


Weight ] 


Fractions of Oxygen 


Y_ = 0.23 
O.o 


Y 0.o= °- 37 


Y n o= °- 50 
U,o 


Y 0 = 0.70 
O,o 


Y„ =0.90 
0,0 












N- Heptane 


8. 35xl0’ 3 


13. 80 x 10" 3 


1 5. 01 x lO* 3 


18. 56xl0’ 3 


22. 34x 10~ 3 


Ethyl 

Alcohol 


8. 58xl0" 3 


11. 1 1 x 1 0* 3 


15. 04xl0* 3 


18.41 xlO" 3 


20. 36xl0" 3 


Benzene 


9.85xl0* 3 


CARBON FORMATION 


Toluene 


7. 70 x 10” 3 


CARBON FORMATION 



on flame shape, and which assumes a discontinuous reaction surface with- 
out proper allowance for dissociation. 

C. Carbon Formation 

For oxygen weight fractions exceeding Yq q = 0. 23, formation of 

carbon residue was obtained during the burning of benzene and toluene. 

Photographs of toluene burning at Y„ = 0. 50 are shown in Fig. 4(b). 

Reference to Fig. 4(b) shows that the carbon formed during combustion at 

Y^ = 0. 50 remained on the filament after all of the liquid fuel had 
u, o 

burned. This carbon residue was inspected and found to be of a dry, 
brittle, consistency. An effort was made to determine the value of Yq^ q 
at which carbon formation is first noticeable. However, no definite con- 
clusions could be drawn since any reasonable change in the oxygen concen- 
tration above that of normal air resulted in carbon formation. The forma- 
tion of carbon in heterogeneous combustion of aromatic fuels has been 
observed by other investigators and merits additional quantitative study. 
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For oxygen weight fraction Y = 0.90, transitory formation of 
carbon occurs during the early stages of combustion of droplets of 
benzene and toluene, but the carbon is consumed by the time combustion 
is completed, i.e. , there are no deposits of the residue remaining on 
the filament at the end of burning. 

Because of carbon formation during the combustion of benzene 
and toluene, it was impossible to determine accurately the burning rates 
of these fuels. The shroudlike formation of carbon caused considerable 
distortion of the drop, and very shortly after ignition the carbon com- 
pletely surrounded and obscured the drop of fuel. 

D. Possibility of Utilizing Diameter Vs. Time Data for Studying the 
Ignition Process 

The experimental data yield not only information on steady 
burning but provide information also on the transitory ignition processes 
prior to steady burning (Fig. 9). Thus a change of burning rate takes 
place during the time that the temperature of the liquid drop is raised 
from its initial ambient temperature to the steady-burning temperature. 
It is possible to study the pre-steady state burning in a manner similar 
to that employed for determining the steady burning rates of drops of 
fuel. However, spark ignition should not be used for the study of 
ignition phenomena. With the electric spark, there is considerable dis- 
tortion of the drop at the time of ignition and, in extreme cases, there 
is some loss of the liquid fuel. 
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IV. COMPARISON BETWEEN THEORY AND EXPERIMENT FOR 
THE EVAPORATION CONSTANT IN OXYGEN-ENRICHED 
NITROGEN-OXYGEN MIXTURES 

The calculated and observed evaporation constants for different 
oxygen weight fractions are contrasted in Table 4. The satisfactory 
agreement between calculated and observed values of K* suggests that 
the physical model upon which the theoretical analysis is based re- 
presents a useful first approximation for predicting the burning rates 
of single droplets of the fuels considered. The success of the theory 
must be ascribed to a nearly complete compensation between obvious 
deficiencies in the analysis: the temperature gradient near the drop 
surface which determines the burning rate, will probably be similar 
to the calculated profile if proper allowance is made for distributed 
chemical reactions, dissociation, and reassociation. 

Further comparison between theory and experiment is shown 
in Figs. 10 and 11 for n-heptane and ethyl alcohol. These plots indicate 
the calculated trend of the evaporation constant, K', (and hence of m^.) 
with increasing oxygen concentration in the oxidizing medium, and also 
the satisfactory agreement between theoretical and experimentally 



determined values. 



TABLE 4. COMPARISON BETWEEN THEORETICAL AND EXPERIMENTAL VALUES OF K’ 
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V. SUGGESTIONS FOR FUTURE WORK 

Additional work on the burning of single droplets of fuel in an 
oxidizing atmosphere should include the following problems: 

1. A detailed investigation of carbon formation during 
heterogeneous burning of aromatic chemicals. 

2. Comparison between theory and experiment for the 
effects of pressure on the burning rates of liquid fuels. 

3. Effect of temperature of the oxidizing atmosphere on 
the burning rate of liquid fuels. 

4. A detailed study of ignition phenomena. 

5. A detailed study of the effects of convection and 
turbulence on the gaseous phase of the fuel. 
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FIGURE 1. SCHEMATIC DIAGRAM OF BURNING FUEL DROP 
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FIGURE 2. 




SCHEMATIC DIAGRAM OF FLAME FRONT 
SURROUNDING A BURNING DROPLET 
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FIGURE 3. PHOTOGRAPH OF APPARATUS. 
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4a. N-Heptane burning 4b. Toluene burning 
at ^,0=0.50. at Y 0 ,o-°-5°- 



FIGURE 4. PHOTOGRAPHS OF BURNING DROPS OF FUEL 
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FIGURE 5. Plot of a* Vs. t for N-Heotane and Yq o -0.50 

(K 1 =15.01 x lCf cm*sec.') 
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FIGURE 6. 





- 27 - 




t (sec) 

FIGURE 7. Plot of d* Vs. t for Toluene and Yq,o =0 * 2 5 

(K*= 7.70 x 10 s cm* sec!) 
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FIGURE 8. Plot of d* Vs. t for Benzene and Y 0>o - 0 * 2 3 

(K’=9.85 x 10* 3 cmtsec.') 
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FIGURE 9. EXAMPLE OF TRANSITORY IGNITION PROCESS. 

(ETHYL ALCOHOL BURNING IN AIR) 



) 
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FIGURE 10. COMPARISON CF THEORETICAL AND EXPERIMENTAL 
VALUES OF K' FOR K-HEPTANE. 
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FIGURE 11 



COMPARISON OF THEORETICAL AND EXPERIMENTAL 
VALUES OF K' FOR ETHYL ALCOHOL. 
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